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ABSTRACT 


The  purpose  of  this  thesis  is  to  find  the  approximate  distri¬ 
bution  of  sample  regression  coefficients  based  on  samples  from  two 
linear  time  series.  Chapter  I  introduces  the  problem.  In  Chapter  II 
the  distribution  of  the  regression  coefficients  with  known  and  unknown 
means,  when  sampling  from  a  bivariate  normal  population  with  zero 
correlation,  are  derived,  using  Daniels's  (1956)  form  of  Geary's 
(1944)  extension  of  Cramer's  theorem.  In  Chapter  III  Daniels's 
method  is  again  used,  together  with  the  saddle -point  approximation 
and  the  method  of  steepest  descents  to  find  an  asymptotic  approximation 
of  the  distributions  of  the  same  coefficients  as  in  Chapter  II  ,  when 


sampling  from  two  linear  time  series. 
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CHAPTER  I 

INTRODUCTION 


In  regression  and  correlation  analysis  of  bivariate  data,  a 
common  error  is  to  infer  the  existence  of  a  direct  relationship  between 
the  two  variables,  whereas,  in  reality,  both  are  related  to  a  third, 
hidden  variable. 

Let  (x,  ,  y  ),  (x  ,  y  ,  (x  ,  y  ),  be  a  sample  of  size  n  from 

11  2  2  n  n 

a  bivariate  population.  Denote  the  sample  regression  coefficient  of 


y  on  x  by 


n 


^  (x  -x)(y  -y) 


x  i=l 


21 


n 


(1.  1) 


E'-, 


x)‘ 


i=l 


and  the  sample  regression  coefficient  with  known  means  by 


n 


21  x. 


i=l 


,yi 


21  n 


(1.  2) 


I  X, 

*- -  i 


i=l 


The  distribution  of  b  ’  w^en  sampling  from  a  bivariate  normal 
population,  was  first  derived  by  K.  Pearson  (1926)  and  Romanovsky 
(1926). 


In  Chapter  II  of  this  thesis  we  derive  the  corresponding 


X 


-  2 


distributions  of  and  b  in  the  case  where  the  population  has 
zero  correlation,  using  Daniels's  (1956)  form  of  Geary's  (1944)  ex¬ 
tension  of  Cramer's  theorem. 

In  Chapter  III,  section  I,  we  derive  the  approximate  distri¬ 
bution  of  b  when  sampling  at  times  t  ,t  ,..jt  ,  from  two  linear 
time  series,  passing  through  the  origin,  with  normal  residuals.  Thus 


we  have 


x.  =  B,  t.  +  e . , 
i  lii 


y.  =  B_  t.  +  f.,  for  i  =  1,  2,  ...  n, 

l  2  i  i 


(1.3) 


where  and  are  constants,  and  ^e.,  f.  ^  are  independent  N(0,  1) 

random  variables.  We  again  use  Daniels's  method  to  find  an  asymptotic 

approximation  of  the  inversion  integral  (2.  1)  involved  in  using  the 

saddlepoint  approximation  and  the  method  of  steepest  descents  (see 

De  Bruijn  (1961)).  In  section  II  the  same  procedure  is  followed  to 

x 

derive  the  approximate  distribution  of  b  when  sampling  at  times 

Ca  -L 

t  ,  t  .  .  .  ,  t  ,  from  the  series 
12  n 


x .  —  a?  +  B  t.  +  e.  , 
l  1  lii 


y.  =  a  +  B  t.  +  f.,  for  i  =  1,2,.  .  ,  n, 
i  2  2  l  l 


(1.4) 


where  a  and  av  are  constants  and  B  .  B„,  and  s  e.,  f.'V  are  as 

12  12  ^  l  i' 

4 

defined  for  (1.  3).  The  corresponding  distributions  of  b  and  b  , 
the  coefficients  of  regression  of  x  on  y,  can  be  found  by  permutating 
the  indices  in  the  distributions  of  b  and  b  ^  ^  respectively. 


' 


-  3  - 


Daniels's  method  is  given  in  Appendix  I. 

The  distribution  of  b  is  tabulated  for  various  exemplary 

21  x  7 

values  of  B  ,  B^,  and  n,  in  Appendix  II. 


f 


-  4  - 


CHAPTER  II 


THE  DISTRIBUTION  OF  THE  SAMPLE  REGRESSION  COEFFICIENT 


If  r  is  a  statistic  of  the  form  c/c  ,  where  c  is  non-negative, 
the  probability  density  function  (p.  d.  f.  )  of  r  is  given  by 


h  (r )  = 


_JL  /  fa-rT,  T) 

2 Tri  /  J)u 


dT  , 


(2.  1) 


u=0 


T  c  +Tc 

where  M(T  ,  T)  =  E  e  ,  is  the  joint  moment-generating  function 

of  c  and  c.  Integration  is  along  the  imaginary  axis  in  the  complex  T 

plane,  or  any  allowable  deformation.  This  is  Daniels's  form  of 

Geary's  extension  of  Cramer's  theorem. 

Let  (x  ,  y  ),  (x  ,  y  (x  ,  y  ),  be  n  observed  pairs  of 

1  l  L-i  Li  n  n 

values  from  a  bivariate  normal  distribution  with  zero  correlation. 

The  joint  distribution  of  x.  ,x  ,  .  .  .,x  ,y,,yv  ■  •  vy  .  is 
J  12  n  1  2  n 


n 


dF  =  (2tt) 


>’nexp{- HZ  xi 


i=l 


n 


1 1,1 } 


i=l 


(2.  2) 


X  dx  .  .  .  dx  dy  .  .  .  dy 
1  n  1  n 


Let 


n 


/  (x.  -x ) (y .  -y ),  and 

‘  —  i  i 


i=l 


(2.  3) 


C  ' 


> 


-  5  - 


where 


n 

cn  =  (x.-x)2  , 

O  u -  i 

i=  1 

n 


x 


= 1  y  x. , 

n  ( —  i 


and 


i=l 


n 


=  i  y . 

n  <• —  1 1 


i=  1 


The  joint  moment-generating  function  of  c  and  c  is 

o 


oo  oo 


n 


M(T  ,  T)  =  (2tt) 
o 


-n 


exp  * 


-HI  v 


00  -00 


i=l 


n 


+ 


n 

V 

2  T  /  (x.-x)2 


n 


2  T 


o 


i=l 


i=l 


L 


i=l 


-  -11 

(x.-x)(y.-y) 


J 


dx,  .  .  .  dx  dy  .  .  .  dy 
1  n  1  n 


We  now  replace  x^,.  .  ,  x^  by  new  variables  w^.  .  ,  w  by  means  of 


n 


1 


n 


an  orthogonal  transformation  such  that 


Wj  =  Vn  x  , 


(2.4) 


and  apply  a  transformation  with  the  same  matrix  to  y  ,  y  ,  which 

1  n 

are  thus  replaced  by  new  variables  z,,.  .  ,  z  ,  such  that 

In 


j  y  . 


We  have  then 


-  6  - 


n 


> 


and 


(x.  -x)  (y.  -y ) 
1  1 


i=l 


i=l 


xiyi 


-  n  x  y  = 


i=2 


w.  z 
l 


Hence 


M(T  ,  T) 
o 


+  z 


2 


+  (1  -  2  T  ) 
o 


i=2 


w.  2  -  2  T 

l 


i=2 


w. 

l 


Z. 

1 


This  is  the  product  of  n  double  integrals. 


Thus 


f 


where 


A  = 


-  7  - 

(n-1) 

,  T)  - 

Ai  2  • 

o 

—  J 

1  -  2T 

-T 

o 

-T 

1 

Hence 


(n-1) 


M(T  ,  T)  =  (1  -  2T  -  T  ) 
o  o 


2  '  2 


Making  the  substitution 


x 


u  -  T  +  b  T  , 
o  21 


we  obtain 


(n-1) 


x 


M(u-b  T,  T)  =  (1  -  2u  +  2  b  T  -  T  ) 


Differentiating  with  respect  to  u,  and  putting  u  =  0,  we  obtain, 


x 


3  M  (u  -  b21  T,  T) 

3  u 


u=0 


(n+1) 


=  (n-1)  (1+2  b21X  T-T2) 


X 


Thus  from  (2.  1 ),  the  p.  d.  f.  of  b  is  given  by 


x  (n-  1 )  , 

h(b2i  >=  i^r1  (1 


_  (n+ 1 ) 

2  b  XT-T2)  2  dT 


(n+1) 


(2.  5) 


(2.  6) 


The  integrand  ^(T)  =  (1  +  2  b^  T-T2) 


2 


has  two  singularities  at 


-  8  - 


the  points 


T=b21%  Vl  +b2i: 


For  any  real  value  of  these  two  points  are  real  and  are  located 

x 

on  either  side  of  the  origin  and  of  the  point  (b  ,  0)  in  the  complex 

-L 

T  plane.  Thus  we  can  deform  the  path  of  integration  from  the  imaginary 
axis  in  the  complex  T  plane  to  the  para.llel  path  through  the  point 


(b  ,  0)  on  the  real  axis.  The  value  of  the  integral  taken  along  this 


new  path  will  be  unaltered  since 

x 


lim 
real  t— >oo! 


<*21 


(0,  -it) 


(0, it) 

^(T)dT  +/  ^(T)dT 

(b21x,it) 


=  0 


(2.  7) 


where  integration  is  along  the  straight  lines  joining  the  points  indicated. 

x 

Thus  letting  T  =  b  +  i 7J  in  (2.  6)  we  have 


oo 


(n+1) 


h<b21X) 


(n-1) 


m 


/  l>b 


21 


+  7  1 


(2.  8) 


Let 


1 


1  + 


V 


(1+b 


x 


21 


] 


-1 


so  that 


d  -q  -  - 


<1+b21X2>2 


2%  (l--g) 


d-e 


: 


' 


-  9  - 


Thus 


(n-1) 

2-rr 


x 


(1  +b21  > 


n 

2 


/* 


s.i  .1 
2  (i- s )  2  a-g , 


(2.  9) 


which  reduces  to 


hfb21S>  = 


r<t) 


-A  r  ) 


(1+b21  2> 


n 

2 


(2.  10) 


If  we  define  the  new  variable 


x 


t  =  ^/n- 1  b 


21  ’ 


it  can  be  seen  that  the  p.  d.  f.  of  t  is  given  by 


P(t  )  = 


r<7> 


V^(n-l)  ) 


n 

t*2  ~2 

(  1  +  —7) 


n-1.  '  n-1 


(2.  11) 


This  is  Student's  't1  distribution  with  n-1  degrees  of  freedom.  If 
we  replace  (2.  3)  by 


n 


Z,T. 

/ - -  1  1 


c  =  /  x.v.  and 

i=l 


(2.  12) 


n 


y  , 

-  /  X.  , 

o  £- — -  1 


i=l 

and  follow  the  same  procedure  without  the  orthogonal  transformation 


(2.  4)  we  see  that  the  p.  d.  f.  of  b  is  given  by 


10  - 


9 


h(b 


21 


(1  +  b 


21 


(n+1) 

2 


(2.  13) 


Thus  the  variable 


t  —  ^ 


b 


21 


has  Student's  1 1'  distribution  with  n  degrees  of  freedom. 
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CHAPTER  III 

THE  APPROXIMATE  DISTRIBUTION  OF  THE  SAMPLE  REGRESSION 
COEFFICIENT  BASED  ON  SAMPLES  FROM  TWO  LINEAR  TIME 

SERIES 


Section  1 

In  this  section  we  follow  the  method  of  the  second  chapter 
to  find  the  approximate  distribution  of 

t  Vi 

b  =  — - 

21  n 


i=l 


where 


x.  =  B  t.  +  e .  , 
l  lii 


(3.  1.  1) 


v=B  t.  +  f for  i  =  1 , 2,  .  .  n  , 
7  i  2  i  i 


where  B}  and  B^  are  constants  and  s^e.,  f.J.  are  independent  N(0,  1) 


random  variables. 


We  consider  the  case 


n 


H  ti2  =  °(n)’ 


(3.  1.  2) 


i=l 


If,  for  example,  we  choose  a  fixed  time  interval  of  length  K  and  take 


12  - 


the  sample  values  at  n  equally  spaced  time  throughout  this  interval, 
we  have 


t.  =  K  —  for  i  =  1 , 2,  .  .  . ,  n. 
l  n 


n 


This  gives  /  t.  2  =  7 —  (n+1)  (2n+l)  =  O(n)  . 
— «  1  on 

i=l 


We  consider  this  case  to  be  more  interesting  than  any  others,  as  it 

would  probably  be  of  more  practical  value. 

The  joint  distribution  of  x,  ,  x_  .  .  .t  x  ,  y.  ,  y y  ,  as 

12  n  1  2  n 

defined  in  (3.  1.  1)  is 


n 

dF  =  (2it)  n  exp  ^  ^  (x.  -  B  ^t  J 2 

i=l 


Let 


+ 


n 


H  (^i  -  Vi’ 


dx  .  .  .  dx  dy  . . .  dy  . 
1  n  1  n 


and 


(3. 1.  3) 


The  joint  moment-generating  function  of  c^  and  c  is 


00  00 


n 


M(T  ,  T)  =  (2tt) 
o 


-n 


J  ■  Wi 

-00  -00  i-1 


.«  •  %  tb  n*b  "  w  "I 


t 


-  13  - 


n 

+  J~  (y.-B  t.)2  -  2T  J  x . 


n 


i=  1 


i=l 


(3.  1.4) 


n 


r-r. 

•T  1  l 


2T 

i=  1 
oo  oo 


1 


'r  dx  .  .  .  dx  dy  .  .  .  dy  , 
1  n  1  n 


=  (2") 


"  J  ■■■  f  exp {’  \ [ 


n 


(1  -2T  ) 
o 


L*.- 

i 


OO  -  00 


i=l 


n  n 

-  2  B  £ _  x.t.  +  y"  y.2 

1  ii  i 

i=l  i=l 


-  2  T 


x.y.  +  (B  2  +  B  2 ) 

ii  1  2 


i=l 


dx_  .  .  .  dx  dy  .  .  .  dy 
1  n  1  n 


Except  for  a  constant  factor,  (3.  1.4)  is  the  product  of  n  double  inte¬ 
grals  of  the  form 
oo  oo 

exp  (-1-  (x,  y)  (*)  -  2  (*)  j  j-  dx  dy  ,  (3.1.  5) 

-  oo  -  oo 


where  A  is  a  symmetric  2x2  constant  matrix,  and  a  and  b  are  con¬ 


stants.  We  denote  the  integrand  of  (3.  1 .  5 )  by 


r  1 


exp  -j  -  -  ^(x,  y) 


I 
J  ’ 


(3.  1.  6) 


where  d(x,y)  can  be  written  in  the  form 


14  - 


^(x,  y)  =  (x,  y) 


[*  ?>  ■  o] 


+  (a 


,b)[— 1 


(*)  " 
b  y 


(a,  b)  A  1  (*) 


(3.  1.  7) 


or 


J^(x,  y)  =  |  (x,  y)  -  (a,  b)  A~  \ 


A 


-  (a,b)  A'1  (b  . 
—  b 


Since  A  is  symmetric,  we  can  write  (3.  1.  8)  as 

(w,  z)  A  (W)  -  (a,  b)  A  1  (f")  , 

—  z  —  b 

where 


(3.  1.  8) 


(w,  z)  =  (x,  y)  -  (a,  b)  A  1  . 


(3.  1.9) 


Clearly  the  Jacobian  of  this  transformation  is  unity,  so 
oo  oo 

/  /  exp  j"-  —  jz4(x,  y)  | dx  dy 


oo  -  oo 


exp  jj-  (a,  b )  A  1  (^)J 


oo  oo 


X 


f  I  exp  -|  -  7  (w,  z)  A  (^^dw  dz 


r 


-  oo  -  oo 


(3.  1.  10) 


2  7T 


I  a! 


exp  ^  -  (a,  b)  A  (fe)  j- 


Using  this  result  to  evaluate  (3.  1.4)  we  have 


, 


15  - 


n 


M(T  ,  T)  =  (2tt)' 


n 


IT 


i=l 


"exp-f-  |  ^  t.2 

i=  1 


'~^r  “pfi  Vi  1  ol 

!a|2 


(3.  1.  11) 


where 


Hence 


A  = 

1-2T 

o 

-T 

-T 

1 

a.  =  B.  t.  ,  and  b  =  B_ 

t.  , 

l  1  i 

l  2 

l 

for  i  =  1,2,... 

i  n. 

i 

A"1  =  (1  -  2T 

-T2)"1 

| 

—  o 

i  t 

1  -2T 


and  we  obtain 


n 


M(T  ,  T) 
o 


(1  - 2T  -  T  2  )  exp 
o 


n 


{ i  i 


■+BZ  > 


i=l 


x  exp  ' 


L 


i=l 


2  (1  -  2T  -  T  ) 
o 


[b}  2+(1-2Tq)B22  +  2T  BiB2fjr 


(3.  1.12) 


' 


. 


9 


-  16  - 


Putting  m  = 


n 


t.  2  and 

l 


i=l 


u  =  T  +  b  T  , 
o  21 


we  obtain 


f  (Bi2+Bz2) 


M(u-b21  T,  T ) 


(l-2u  +  2b21  T-T2) 


n 

2 


r 


m 


x  ^tTT-Zu.Zb  T-T2)  IV  +  (1-2u+2hiT)  b2  +2TB1B2 

21 

(3.  1.  13) 


r 


Differentiating  with  respect  to  u  we  have 

d  M  (u-b  T,  T) 


[- 


n 


M (u-b 2i  T,  T)  |  (f_2u+2b21  T-T2) 


+ 


m 


r 


-  B  (1  - 2u+2b  T-T  ) 


(l-2u+2b21  T-T2)2  |  ~2  "21 


+ 


[V+(1 


-2u+2b  T)  B  "+2T  B.B.ll. 


21  7  2 


T 


2lr 


(3.  1.  14) 


and  putting  u  =  0  we  have 


17  - 


(u-b  T,  T) 


u=0 


[' 


<1+2b21  T'T  > 


n  (  n(l  +  2b2J  T-T")  +  m(B2T+B1)a 
2.2  +Z 


x  exp 


B,2  +  BZ»  +  2TB  <b  B  +B  j 


{?[“ - i 


+  zb21  T‘T 


J.  (3.1.1 


5) 


Thus  from  (2.  1 )  the  p.  d.  f.  of  b  is  given  by 


f  <bi2+b22) 


h(b,,) 


21  2-tri 


;[n(l  +  2b21  T-T2)  +  m(B7T+B1)2J 


x 


[ 


1  +  2b21  T-T 


] 


f+2 


u 

x  exp-^Y 


Bl2+B22+2TB2(b21B2+Bl)ll 


1  +  2b  T-T 

n  L*  1 


JJ 


dT  , 


(3.  1.16) 


where  integration  is  along  the  imaginary  axis  in  the  complex  T  plane, 
or  along  any  allowable  deformation. 

In  the  case  (3.  1.  2),  with  a  possible  relative  error  which  is 


-1 


0(n~  )  we  may  write  m  =  Rn,  where  R  is  a  constant.  In  (3.  16)  where 


m  occurs  it  is  multiplied  by  homogeneous  terms  of  the  second  degree 

in  B  and  B  and  we  may,  without  loss  of  generality,  take 

1  2 


m  =  n  , 


. 


. 


18  - 


and  write  (3.  1.  16)  in  the  form 


n 


"7Bi 

,  n  x _ n  2  1 

^21  2  ui  6 


C1+2b21T“T2  +  (B2T+Bi)2] 


n 


[ 


1  +  2b21  T-T 


1 


+  2 


x  exp 


{■ 


nr  Wl,, 
l+2b21T-T2  ,K  C 


Let 


z  = 


T  '  b21  .  b2lVBl 

-  ,  A  =  - 


Vl+b 


21 


^1+V 


p2  =  2  +  A2  ,  and  k  =  B^A  . 


Then  we  have 


(1  +  2  b21  T-T2)  -  (1  +  b212)  -  (b21-T)2 


-  (1  +  b212)  (1  -  z2)  , 


and 


(B. 


T  +  B  i ) 2  - 


=  (1  +  b212)  (B2z  +  A)* 


Thus 


(1  +  2b  T-T2)  +  (B2T  +  B^)' 


=  (1  +  b2J  2)  |*1  +  A2  +  2kz  +  (B22-1)z2]  , 


(3.  1.  17) 


(3.  1.18) 


and 


» 
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(B2  T  +  B  1 ) 2 
~('l+2b  T-T2) 


B  2z 2  +  2kz  +  A2 
(1  -  z2) 


p2  z2  +  2kz  2 

(1-z2)  +A 


Substituting  in  (3.  1.  17)  we  have 


-|(B  2-A2)  J&l) 

h(b21  >  =  ATe  <1+b212> 


[l+A2  +  2kz  +  (B  2-l)z2] 
(1  -  z2)2 


x  exp  * 


n 

2 


2  2 


2kz  +  p  z 

2 


1  -  Z 


■/^n(l  -  z  2) 


1 


r  dz  , 


I 


(3.  1.  19) 


where  integration  is  along  the  path 


Re  (z )  =  - 


21 


vi+b2i2 


(3.  1.  20) 


or  any  allowable  deformation.  We  evaluate  this  integral  using  the 
following  saddlepoint  approximation  and  method  of  steepest  descents, 


/ 


f(x)etg(x)dx 


=  V2^  “  [t  |g1 1  (x)|  J  ■  f(x) 


tg(x) 


X 


I  1  +  0(t  1 )}  (t  — >oo )  ,  (3.  1.  21) 
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where 


A 

x=x 


(3.  1.  22) 


A  dg 

x  is  the  solution  of  — =  0  , 

dx 

the  path  of  integration  in  the  neighbourhood  of  the  saddlepoint  at  x 
is  the  path  of  steepest  descent,  which  is  directed  along  the  axis  of 
the  saddlepoint,  that  is,  along  the  straight  line 

(x-x)2  g  (&)  real  and  <  0  ,  (3.1.23) 


and  a  is  the  complex  number  with  modulus  1  whose  argument  corres¬ 
ponds  to  the  direction  of  the  axis  (see,  for  example,  De  Bruijn  (1961), 
Chapter  5 ). 


The  integral  in  (3.  1.  19)  is  in  the  form  of  the  left  member 
of  (3.  1.21)  with  x  replaced  by  z,  where 


,  ,  2kz  +  p2z2  0  /i  2\ 

g(z)  =  — - - S -  -  An(l  -  z  ), 

J-  —  z 

[l  +  A2  +  2kz  +  (B22-1)z2J 


f  (z)  = 


(1  -  z2)2 


(3.  1.  24) 


n 


and  t  =  —  . 


The  only  singularities  of  the  integrand  of  (3.  1.19)  are  at 

z  =  +  1 

and,  from  (3.  1.  20),  for  any  finite  value  of  b^,  the  path  of  integration 


' 


' 
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in  the  complex  z  plane  is  parallel  to  the  imaginary  axis  and  cuts 


the  real  axis  at  a  point  in  the  interval  (-1,  1).  We  shall  show  that 


there  is  always  a  unique  point  z  £  (-1,1),  such  that 


dg.  I 

dz  I 

*  A 

z  —  z 


From  (3.  1.  24)  we  have 

df  "JT-P?  0+  (p2+l)z  +  kz2  -  z3]  (3.1.25) 


and 


(irby  Cp2+1  +  2kz  ■ 3z^ 


4z 

(1  -z2) 


dz 


Thus,  from  (3.  1.  22), 


gll*®)  (1,SV 


0*' 


+  1  +  2kz  -  3z2  J 


(3.  1.  26) 


Let 


h(z )  =  k  +  (p  2+l  )z  +  kz  2  -  z  . 


From  (3.  1.  25)  and  (3.  1.  26)  we  see  that  for  real  z,  if  z2  /  1, 


h(z) 


=  0 


0 


(3.  1.  27) 


A 

z  =  z 


and 


dh 


o  (z)  has  the  same  sign  as  — 
&llv  ’  &  d; 


A 

z  =  z 


(3.  1.  28) 


■ 


■ 
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dh  . 

Also  h(z)  is  a  continuous  function  of  z  with  three  roots  and  — —  is 

dz 

a  continuous  function  of  z  with  two  roots.  We  denote  the  real  parts 
of  the  roots  of  h(z)  by  A  <  B  <  C. 


From  the  definitions  (3.  1.18)  we  note 


p2  >  f  2k  | 


(a)  Suppose  p2  ^  +  2k  . 

Then  p 2  >  |  2k  | 


and  we  have 


h(-l )  =  2k  -  p2  <  0  , 
h(l )  =  2k  +  p2  >  0  , 


h(oo)  =  -  oo  and 


h(-oo)  -  oo  . 


Thus  h (z )  has  three  real  roots  where 


A  £  (-oo,  -1),  B  £  (-1,  1)  ,  and  C  G  (1,  oo). 


Let  z  =  B . 


(b)  Suppose  p2  =  21c  4  0,  which  implies  that  k  >  0.  Then 


h(-l )  =  0  ,  h(l)  =  4k  >  0  , 


h (- oo)=  oo  ,  h(oo)=  -oo  ,  and 


dh 

dz 


=  -2  <  0  . 


Thus  h(z )  has  three  real  roots  where 

A  =  - 1 ,  B  £(-1,1),  and  C  G(l,oo). 


Let  z  =  B . 


. 
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(c) 


Suppose  p2  =  -2k  1  0,  which  implies  that  k  <  0.  Then 
h(-  1 )  =  4k  <  0  ,  h(l)  =  0  , 

h(-oo)  =  oo  ,  h(oo)  =  - oo  ,  and 


dh 

dz 

z  =  1 


[ 


2k(z  -  1 )  +  (1  -3z  ) 


l 


z  =  1 


=  -2  <  0  . 


Thus  h(z)  has  three  real  roots  where 

A  £  ( -  oo ,  -1)  ,  B  £  (-1,  1),  and  C  =  1. 
Let  'z,  =  B. 

(d)  Suppose  k  =  0  .  Then 

h(z)  =  z(p2+l  -  z2)  , 


so 


h(-l)  =  -p2  <  0  ,  h(0)  =  0,  and 

h(l)  =  p2  >  0  . 

Thus  h(z)  has  three  real  roots  where 


AS  (-oo,  - 1^|  ,  B  =  .0  ,  and  C  E  j^l,  oo) 


Let  z  =  B. 

In  all  four  cases  (a),  (b),  (c),  (d),  — —  must  have  a  root  in 

dz 

each  of  the  intervals  (A,  B)  and  (B,C),  so 


dh 

dz 


t  0  > 


A 

z  =  z 


and  since  h(z)  <  0  for  z  £  (A,  B)  ,  and  h(z)  >  0  for  z  £  (B,  C),  we 


have 


» 
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dh 

dz 


A 

z  =  z 


>  0  . 


Thus  from  (3.  1.  27)  and  (3.  1.  28)  we  see  that  for  any  finite  values 
of  p2  and  k  there  exists  a  unique  real  number  z  such  that 


Z  £  (-1,  1)  ,  ^ 
dz 


=  0  ,  and  g  1  (£)  >  0. 


A 
Z  =  Z 


F r om  (3.  1 .  23 )  we  see  that,  since  z  is  real,  and  g  ^  (z )  >  0, 
the  axis  of  the  saddlepoint  is  the  straight  line 

Re  (z )  =  z  . 

Thus,  in  the  neighbourhood  of  the  saddlepoint  we  can  deform  the 
path  of  integration  of  (3.  1.  19)  to  pass  along  the  axis  of  the  saddle¬ 
point,  and  use  the  saddlepoint  approximation  (3.  1.21)  with  (3.  1.  26), 
to  obtain 


h(b 


21 


(1  +  b 


21 


(n+1) 

2 


[  1  +  A2  +  2kz  +  (B22-1)z2‘] 

x - T 

n  ,  1 
yd  1  T  ? 

(1-z2)  J”p2  +  l+2kz-3z2j 

x  |  1  +  O (n  *)  y  (n  ->  oo).  (3.  1.  29) 

Putting  =  B^  =  0  in  definitions  (3.  1.18)  we  have  A  =  p2  =  k  =  0, 
hence  from  case  (d)  we  have  z  =  0.  Substituting  these  values  in 


exp 


2  k"z  +  p  2 

l2 

.  o  -  J 

ll 


J 


25  - 


(3.  1.29)  we  obtain 


h(b. 


i)  =  Jt  <1+b2d> 


(n+1) 


2.  2 


x  |  1  +  0(n  1 )  } 


(n  ->  oo) 


(3.  1.30) 


As  expected  (3.  1.  30  )  is  consistent  with  (2.  13)  which  can  be  obtained 
exactly  from  (3.  1 .  3  0)  by  renormalising. 


Section  II 


In  this  section  we  follow  the  procedure  of  Section  I  to  find 

the  approximate  distribution  of 

n 

\  (x.-x)  (y.-y) 

4 _  i  i 

x  i=  1 


21  n  , 

Y_  ^  -  x)2 

i=i 


where 


x.  =  a  +  B  t.  +  e.  , 
i  1  1  i  i 


v  =  +  B  t  +  f .  ,  for  i  =  1 , 2,  .  .  . ,  n  , 

}  i  2  2  i  i 


(3.2.  1) 


where  are  constants,  and  are  independent 

N(0,  1)  random  variables. 


The  joint  distribution  of  x^ ,  x^,  .  .  .,  x^,  y^  ,  y^,  .  •  • ,  y^, 


as 


defined  in  (3.  2.  1)  is 


. 


i  .  *  '4-1 1 
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dF 


-  B  t  )2 
1  1 


+ 


n 

H 

i=l 


(y,  - 


a. 


B  t. ) 

2  l 


l 

I 


dx  . . .  dx  dy  . .  .  dy  . 
1  n  1  n 


Let 


and 


n 


i=l 


The  joint  moment-generating  function  of  c  and  c  is 


(3.  2.  2) 


M(T  ,  T)  =  (2tt) 
o 


a,  -  B  t.)2  +  )  (y. 

1  1  i  L —  i 

i=l 


-  B  t. ) 2 

2  l 


(3.  2.  3) 


x  dx  .  .  .  dx  dy  .  .  .  dy  . 
1  n  1  n 
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We  replace  x^, .  .  x^  by  w^, .  .  . ,  w^  by  an  orthogonal  transformation 


w. 


w 


n 


B 


x 


n 


where  B  is  an  n  x  n  orthogonal  matrix  such  that 


w 


!  =  Vn 


x 


Thus 


B  = 


J_  _1 

V  n  V 


n 


_1_ 
V  n 


C 


(3.  2.  4) 


where  C  is  an  (n-  1 )  x  n  matrix.  We  apply  the  same  transformation 


to  y  ,  .  . .  y  to  obtain  %  .  z  , 

1  n  1  .  n 


where 


n 


B 


Y1 


y 


n 


and  thus  from  the  definition  of  B  we  have 


z  =  Vn  y 


We  have  then 


n 

V  , 
/_  x 


n 


n 


=  Z  W  »  , 


1  —  1 
i=l  i=l 


n 

V 

^1  '  1 
i=l  i- 1 


n./. 


(3.  2.  5) 
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n 

L 


(x.-x)‘ 


i=l 


n 


L 


i=2 


\v. 

1 


n  n 

XL  <Vx)  (y.-y) 

i=l 


Lw,  z .  j 

i  l 


n 

7  X  t 
/—  x- 


1  1 


-  (t. 


i=l 


t  ) 
n 


x . 


n 


=  (t  ...  t  )  B 
1  n  — 


-1 


w. 


w 


n 


n 


r — 

and  /  y.  t.  =  (t  • 
* —  ii  I 


t  )  B 

n  — 


i=l 


n 


Substituting  (3.  2.  5)  in  (3.  2.  3)  we  obtain 


M(T  ,  T)  -  (2tt) 

o 


-n 


oo  oo 

I  ^  I 

-00  -OO 


_  n  n 

I  '[V  2  Y 

exp  J  w  2  +  zi 

L  Li.l 


n 


-  2T 


V" 

/  w 
o  f- —  i 


i=2 


n 

2T  /  w.z. 
i=2 


-  2  Vn  -  2  Vn 


2Bl<‘l  •••  ‘n’5. 


-1 


?i 

-  2B2(tl  ... 

-1 

t  )B 
n  — 

zi 

w 

n 

z 

n 

n 


+  n («j2  +  ^2)  +  Y-  S'  <B12+B22)  +  2n‘  (0flBl+“2B2' 


l 


L—  1 


J 

(3.  2.  6) 


x 


dw.  .  .  .  dw  dz,  ...  dz 
1  n  i 


n 
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n 


-  i  V" 

where  t  =  —  /  t 

n  1 — ’  i 

i=l 


Since  B  is  orthogonal  we  have 


B_1  =  B'  = 


1 

1  c 
> 

_1_ 

•y  n 

C' 

J_ 

1  fl 
> 

thus 


-1 


(t,  .  .  .  t  )  B  ‘  =  (Vn  t  I  (t 
1  n  —  '  i 

so  (3.  2.  6)  can  be  written  as 


•  t  )  C')  , 
n  — 


M(TQ,  T)  -  (2 tt ) ~ n  exp  ■£-  |  [n(»1 


+  *2  ) 


(3.  2.  7) 


(3.  2.  8) 


+  (B/  +  B22)  ^  t.2  +  2nt(or1B1  +  W]| 


i=l 


•J 


oo  oo 


X 


f  ...  j  ex p  \  [wp  +  ^2  -  2o'1  VnWj  -  2a  z  yn 


-OO  -OO 


n 


-  2  B  -/n  t  -  2  B2  /n  t  +  (1-2TJ  wi' 


i=2 


+ 


n  n 

2_  z.2  -  2  T  Z_  w.z.  -  2B  j  (t[  ...  tn)  C' 

i=2 


w. 


w 


n 
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-  2B  (t  .  .  .  .  t  )  C’ 
2  1  n  — 


'2 


n 


ll 


dw  .  .  .  dw  dz  .  .  .  dz 
i  n  1  n 


(3.  2.  9) 


Integrating  with  respect  to  Wj  and  z  we  obtain 


M(T  ,  T)  -  (2TT)“(n"1}  e 
o 


n 

1=1 


oo  oo 

-oo  -oo 


exp 


(1 


n  n  n 

•2T  )  )  w.2+  /  z . 2 - 2T  /  w.z. 

O  / -  1  L -  1  -  l  l 


i=2 


i=2 


i=2 


2B  (t  ...  t  )  C' 
11  n  — 


w. 


w 


n 


-  2B2(tl  •• 


t  )  c' 

n  — 


n 


1} 


X  dw.  .  .  .  dw  dz„  ...  dz 
2  n  2  n 


(3.  2.  10) 


Except  for  a  constant  factor  this  is  the  product  of  n-1  double  integrals 
of  the  form  (3.  1.  5).  Thus  from  (3.  1.  10)  , 


n 


1 


M(T  ,  T) 
o 


i<V+V>  b 

i=  l 


-  t) 


=  e 


(n-1) 


n 


A 


exp 


Y  \  (a.b.)  A 

C~~~  C  11  — 


i=2 


a. 

i 

b 

i 


}■ 


(3.  2.  11 ) 


„ 
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where 


A  = 

1-2  T 

o 

-T 

2  ‘  • 

an)  =  Bl(tl  ' 

.  t  )  C' 

n  — 

2  *  ‘  • 

V  =  B2(tl  ' 

.  t  )  C' 

n  — 

Thus 


n 


1 


M(T  ,  T)  =  e 
o 


,(B  2+B  2)  /I  (t.- 
2  1  2  .  ,  1 

1 


t)' 


n 


,  Y*  (a. 2  +(1  -2T  )b.2  +  2Tab.) 

(n-1 )  r  ^—  i  o  1  11 

I  i  i=; 


2 .  2  I  1  i=2 

x  (1  -2T  -T2)  exp  % 

o 


l 


2  (1-2T  - T  ) 

o 


(3.  2.  12) 


where 


n 


n 


1  V 


B 


^  zL  ai2  =  b*2  Z_  V 


i  =  2 


i=2 


B1B2 


n 

V" 

/  a.  b 


i  i 


i=2 


(t  j  •  • 


t  )  C'  c 

n  —  — 


n 


(3.  2.  13) 


Since  B  is  orthogonal  we  have 


4 
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V 

L— 


i=  1 


t  )  B'  B 
n  —  — 


t 

n 


(3.  2.  14) 


thus  from  (3.  2.  8) 


n 


i_  t  |(t  ...  t  )  o 

'  i  In  — 

i=  1 

Vn  t 

c 

5  | 

t 

n 

(3.  2.  15) 


=  n  t2  +  (tj 


t  )  C'  C 
n  —  — 


t 

n 


thus  from  (3.  2.  13) 


V  2 

Z_  ai 


i=2 


n 


i=2 


B1B2 


n 


i=2 


n 


\ 

Z_  (t, 


■t) 


i=l 


(3.  2.  16) 


Substituting  (3.  2.  16)  in  (3.  2.  12)  we  obtain 


M(T  ,  T) 
o 


'2(bi2+b22)  21  (y*)2 


(3.  2.  17) 


. 
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n 


(n-1) 


x  (1 -2T  -T  ) 
o 


2  '  2 


exp 


r  21  (yt  )■ 

■U=1 


^2(l-2To-T^)  lLBl2+(1-2To)B22+2TB1B2]| 


We  see  that  (3.  2.  17)  can  be  obtained  from  (3.  1.12)  by  replacing 


n 

V 


n 


/  t.  and  (1  -  2T  -  T  )  by 
1 —  1  o 

i=  1 


n 


I  V 


(n-1) 

—  ,  ,2 

t)  and  (1-2T^-T  )  respectively. 


i=l 


Thus,  if  we  follow  the  same  procedure  as  in  Section  I  from  (3.  1.  12) 


onwar  d  s ,  with 


m-1  = 


n 

S". 

(t.-t  ‘ 

i=l 


and  assume 


m  =  O(n)  , 


we  obtain  the  expression  corresponding  to  (3.  1.29  )  for  the  p.  d.  f. 
x 


of  b 


21  ’ 


h(b21X> 


V  2tt 


<iliii(Bl2+B22)  x  -f 

c  v  L  n  ,  i  x2\  L 

(1+b21  > 


(3.  2.  18) 


- 


I 
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[l  +  A2  +  2kz  +  (B  2-l  )z2^] 


n-I 

n  A  2  \  2 
(1-z  ) 


+  1 


£p2-rl  +  2kz-3z2J 


1_ 

2 


exp 


2kz  +  p 

(1  -  z  ) 


2 


X 


x  |l  +  0(n  )|  (n— >  oo)  , 

where  A,  k,  p2  are  defined  by  (3.  1.  18)  with  b  replaced  by  b 

c-j  -L  L  1 

A  Fortran  program  was  written  to  compile  re -normalised 

tables  of  the  approximate  distribution  of  b  given  in  (3.  1.  29)-  It 

x 

can  be  seen  from  (3.  2.  18)  that  the  distribution  of  b  for  n  degrees 


of 


freedom  is  the  same  as  that  of  b  for  n- 1  degrees  of  freedom. 


Also  from  the  distribution  of  b  (3.  1.  29),  definitions  (3.  1.  18),  and 


the  fact  that  z  is  a  root  of  (3.  1.  25)  it  can  be  seen  that  the  following 

*2 


pairs  of  values  for  B^,  B9  and  n  give  identical  distributions. 


(i) 

Bi 

=  V 

B2  =  b2' 

n,  and  B^  = 

-v 

B2  =  A’  n; 

(ii) 

Bi 

=  0, 

B^  =  b,  n, 

and  B  =  0, 

ii 

(N3 

PQ 

-b,  n; 

(iii) 

Bi 

=  b, 

B  =  0,  n 

2 

and  B  =  -b, 

B2  = 

0,  n; 

and  that  the  distribution  is  symmetric  if  either  B^  =  0  or  B^  =  0. 
Also  the  values  B  =  ,  B^b^n,  give  a  distribution  which  is  a 

reflection  in  the  axis  of  b  of  the  distribution  given  by  BJ  =  -b, 


B2  =  b2’n' 


The  tables  in  Appendix  II  show  the  5  distinct  distributions  of 
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^21  where  the  values  for  B  ^  and  are  chosen  from  the  integers  -1, 

0,  1,  for  each  of  the  values  15,  30,  45,  for  n.  It  can  be  seen  that  in 

all  the  cases  a  measure  of  central  tendency  is  the  value 

n 


B 


1 


1 +B 1  ^ 


in  the  case  (3.  1.  2),  and  this  value  is  the  mean  in  the  cases  where  the 
distribution  is  symmetric.  In  each  case  the  variance  of  b  decreases 
approximately  as  For  each  value  of  n  the  variance  for  the  case 
B^  =  0,  B^  -  1  is  approximately  twice  that  for  the  case  B^  =  0,  B^  =  0, 
which  latter  is  approximately  twice  that  for  the  remaining  cases  B^  =  1, 


B2  =  0;  B1  =  1,  B2  -  1;  B^=  -1,  B2  -  1. 


* 
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APPENDIX  I 


DANIELS'S  FORM  OF  GEARY'S  (1944)  EXTENSION  OF  CRAMER' S 

THEOREM 


We  are  interested  in  the  distribution  of  statistics  of  the 
c 

form  r  -  ,  where  c  is  non-negative.  If  c  ,  c  have  a  joint  p.  d.  f. 

o 

f(c  ,c),  the  density  for  r  is 
o 


oo 


h(r)  - 


c  f(c  ,  r  c  )  d  c  . 
o  o  o  o 


o 


Let 


M(T  ,  T)  =  E  e 
o 


T  c  +Tc 
o  o 


be  the  joint  moment-generating  function  for  c  and  c.  We  are  con¬ 
cerned  only  with  cases  where  M(T  ,  T)  exists  in  strips  of  non-zero 

o 

width  containing  the  imaginary  axes  in  the  T^  and  T  planes.  The 
usual  Fourier  inversion  formula  is  most  conveniently  written  as 


f  (c  ,  c)  =  —  V  -2 
o  (2m) 


-T  c  -Tc 

M(T  ,  T)  e  °  °  dTQdT  , 


where  integration  is  along  the  imaginary  axes  of  T^  and  T,  or  any 


allowable  deformations  of  these  paths.  In  particular, 
1 


f(c  ,  rc  )  = 
o 


rp  ~(T  +r  T  )c 

,  T  /  /  M(T  ,  T)  e  °  °  d  T  dT 

o7  (2m)2  JJ  o’  ’  o 


rr  "u  c0 

- — r  //  M(U  -  r  T,  T )  e  dudT, 

tri)  ov 


(2-rri ) 


■ 


■ 
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where  the  integration  of  u  =  T  +  rT  is  taken  over  a  similar  path 
in  the  u  plane.  Inversion  of  the  transform  with  respect  to  u  gives 
oo 


r 


o 


u  c 


f(c 


C  ) 
o 


o 


d  c 

o 


1 

Z  rri 


M(u  -  r  T, 


T)  dT  , 


so  that,  when  differentiation  is  permissible, 


oo 


f  c  f(c  , 

J  o  o 


u  c 

r  c  )  e  d  c 
o  o 


—  /■ 


1  /  £M(u  -  r  T,  T) 


dT 


and 


h(r)  - 


i,  f 

.  TTl  J 


3M(u  -  rT,  T) 

3  u 


u-0 


dT  . 
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APPENDIX  II 


EXAMPLE  TABLES  OF  THE  APPROXIMATE  DISTRIBUTION  OF  THE 

SAMPLE  REGRESSION  COEFFICIENT 


The  following  tables  provide  examples  of  the  renormalized 
density  and  distribution  functions  of  the  approximate  distribution  of 
b^.  Thus  if  h(b^^)  is  defined  by  equation  (3.  1.29),  the  columns  in 
the  table  give  values  of  h(b  )/NORM,  H(b  )/NORM  and  z  for  various 

L  1  l  1 

values  of  b  centred  at  B  B  /(I +B  2 ),  where 

Lj  J.  1  L  -L 


oo 


oo 


NORM  = 


h(b2i)db21,  H(b  )  = 


h  (x)dx 


-  oo 


-  oo 


and  z  is  defined  in  section  3.  1.  The  mean  and  variance  of  the  re¬ 
normalized  distribution  and  the  value  NORM  of  the  normalizing  constant 
are  given  at  the  bottom  of  each'table. 


. 
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B1  =  0.  0  B  =  1.  0  N  =  15 


b21 

h(t>2i  )/N°rm 

H(b  )/NORM 

z 

-2,200 

*“.  00004 

.00001 

.39298 

-2.100 

.00007 

.0000] 

.39096 

-2.000 

.00013 

.00002 

. 38866 

-1 .900 

.00022 

.00004 

.38601 

-1.800 

. 00039 

.00007 

.38294 

-1.700 

.00070 

.00012 

.37937 

-1.600 

.00125 

.00022 

.37519 

-1  .500 

.00225 

.00039 

. 37027 

-1 .400 

.00406 

.00070 

.36445 

-1.300 

.00732 

♦  0  0 1 26 

.35751 

-1 .200 

.01314 

.00226 

.34920 

-1.100 

. 02336 

.00405 

.33919 

-1  .000 

. 04094 

.00721 

.32711 

-.900 

.07025 

.01268 

.31248 

I 

9 

o 

o 

.11719 

.02192 

.29478 

-.700 

. 18862 

. 03703 

.27343 

-.600 

. 29057 

. 06079 

.24786 

-.500 

,42513 

.09638 

.21759 

-.400 

.58648 

. 1 4683 

.18232 

-.300 

.75817 

.21407 

.14214 

-.200 

.91 398 

.29787 

.  09758 

-.100 

1 .02386 

.39514 

. 04969 

0.000 

1 .06359 

.50000 

0.00000 

.100 

1 . 02386 

.60486 

-.04969 

.200 

.^1398 

.70213 

-.09758 

.300 

.75817 

.78593 

-.14214 

.400 

. 58648 

. 85317 

-. 18232 

.500 

.42513 

.90362 

-.21759 

.600 

. 29057 

.93921 

-.24786 

.700 

. 1 8862 

.96297 

-.27343 

.300 

.11719 

.97808 

-.29478 

.900 

.07025 

.98732 

-.31248 

1.000 

. 04094 

.99279 

-.3271 1 

1.100 

. 02336 

.99595 

-.33919 

1 .200 

.01314 

.99774 

- . 34920 

1 .300 

.00732 

.99874 

-.35751 

1.400 

.00406 

.99930 

-.36445 

1.500 

. 00225 

. 99961 

-.37027 

1 .600 

.00125 

.99978 

-.37519 

1.700 

.00070 

.99988 

-.37937 

1 .800 

,00039 

. 99993 

-.38294 

1 .900 

.00022 

.99996 

-.38601 

2.000 

.00013 

.99998 

- . 38866 

2.100 

.00007 

.99999 

-.39096 

2.200 

.00004 

.99999 

-.39298 

MEAN  = 

0.00000  VARIANCE  =  .  15327 

NORM  =  1. 

02723 
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B  =  0.  0  B2  =  0.  0  N  =  15 


b21 

h(b2i)/N°RM 

H(t>  )/N0RM 

z 

-2.200 

.00000 

.00000 

0.00000 

-2.100 

.00000 

.00000 

0.00000 

-2.000 

.00000 

.00000 

0.00000 

-1.900 

.00001 

.00000 

0.00000 

-1 .BOO 

.00001 

.00000 

0.00000 

-1  .700 

.00003 

.00000 

0.00000 

-1 .600 

.00006 

.00001 

0.00000 

-1.500 

.00012 

.00002 

0 .00000 

-1 .400 

.00026 

.00004 

0.00000 

-1 .300 

.00055 

.00007 

0.00000 

I 

• 

• 

ro 

O 

o 

.00121 

.00016 

0.00000 

-1.100 

.  00267 

.00035 

0.00000 

o 

o 

c 

. 

1 — 1 

1 

.00594 

.00076 

0.00000 

-.900 

.01319 

.00168 

0.00000 

-.800 

.02904 

.  00372 

0.00000 

-.700 

.06255 

.00815 

0 .00000 

-.600 

.12984 

.01749 

0 .00000 

-.500 

.  25494 

.03628 

0.00000 

-.400 

♦  4  6  3  5 1 

.07162 

0.00000 

-.300 

.76262 

.13242 

0.00000 

-.200 

1  .  11034 

.22602 

0.00000 

o 

o 

r— 1 

• 

1 

1.40330 

. 35245 

0.00000 

0 .000 

1 .51958 

.50000 

0 ,00000 

.  100 

1  .40330 

.64755 

0 .00000 

.200 

1 .11034 

.77398 

0 .00000 

.300 

.76262 

.86758 

0 .00000 

.400 

.46351 

.92838 

0.00000 

.500 

.25494 

.96372 

0 .00000 

.600 

.12984 

.98251 

0 .00000 

.700 

.06255 

.99185 

0.00000 

.80  0 

. 02904 

.99628 

0 .00000 

.900 

.01319 

.99832 

0.00000 

1.000 

.00594 

.99924 

0.00000 

1.100 

. 00267 

.99965 

0 .00000 

1.200 

.00121 

.99984 

0.00000 

1.300 

.00055 

.99993 

0 .00000 

o 

o 

• 

.00026 

.99996 

0.00000 

1.500 

.00012 

.99998 

0.00000 

1 .600 

.00006 

.99999 

0.00000 

1.700 

.00003 

1  .  00000 

0 .00000 

1 .800 

.00001 

1 .00000 

0.00000 

1.900 

.00001 

1 .00000 

0.00000 

2.000 

.00000 

1.00000 

0 .00000 

2.100 

.00000 

1 .00000 

0.00000 

2.200 

.00000 

1 .00000 

0.00000 

mean= 

.00000  VARIANCE^  .07692 

NOR  m  = 

1 .01679 


• 

• 

. 
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B  j  =  1.0  B2  =  0.  0  N  =  15 


b21 

h(b  )/NORM 

H(b  )/NORM 

z 

-1 .500 

.00000 

.00000 

0.00000 

/ 

i 

i — » 

• 

o 

o 

.00000 

.00000 

0.00000 

-1 .300 

.00001 

.00000 

0.00000 

-1 .200 

.00002 

.00000 

0 .00000 

”1.100 

.00005 

.  000  0  0 

0 .00000 

”1 .000 

.00017 

.00002 

0.00000 

-.900 

.00058 

.00005 

0.00000 

-.800 

.00198 

.00017 

0.00000 

-.700 

.00689 

.00057 

0 .00000 

-.600 

.  02358 

.00198 

0.00000 

-.500 

.07662 

.00662 

0.00000 

-.400 

.22567 

.02086 

0.00000 

-.300 

.57070 

.05912 

0.00000 

-.200 

1 .16993 

.14473 

0.00000 

-.100 

1 . 94512 

.29623 

0.00000 

0 .000 

2.15736 

.50000 

0.00000 

.100 

1 .84512 

.70377 

0.00000 

.200 

1 .16993 

.85527 

0 .00000 

.300 

.57070 

. 94088 

0.00000 

.400 

.22567 

.97914 

0 .00000 

.500 

. 07662 

.99338 

0.00000 

.600 

. 02358 

.99802 

0.00000 

.700 

.00689 

.99943 

0.00000 

.800 

.00198 

.99983 

0.00000 

.900 

.00058 

.99995 

0 .00000 

1  .000 

.00017 

.99998 

0.00000 

1.10  0 

.00005 

1.00000 

0.00000 

1.200 

.00002 

1.00000 

0 .00000 

1.300 

.00001 

1.00000 

0.00000 

1 .400 

.00000 

1  .  00000 

0.00000 

1.500 

.00000 

1.00000 

0 .00000 

mean= 

.  00000  VAR  I ANCE=  .  03715 

NOPM  = 

1 .01286 
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b.  -  : 

1.0  B  = 

1.  0 

N  =  15 

l 

2 

b2) 

h(b2l)/N0RM 

H(b  )/NORM 

z 

1  .  150 

.00000 

.00000 

.  0491  5 

1.050 

.00000 

.00000 

.01724 

-.950 

.00000 

.00000 

-.01812 

-.350 

.00001 

.00000 

-.05705 

-.750 

.00002 

.00000 

-.09949 

-.650 

.00005 

.00001 

- . 14510 

-.550 

.00014 

.00001 

-.19309 

-.45  0 

.00041 

.00004 

-.24212 

-.350 

.00128 

.00012 

-.29019 

-.250 

.00419 

.00037 

-.33471 

-.150 

.01419 

,00121 

-.37293 

-.050 

. 04750 

.00406 

-.40273 

.050 

.14787 

.01319 

-.42349 

.150 

. 39832 

.03926 

-.43625 

.250 

. 86894 

.10125 

-.44307 

.350 

1 . 46669 

.21803 

- . 4461 4 

.450 

1 .88258 

.38773 

-.44713 

.550 

1 . 85500 

.57761 

-.44715 

.650 

1 .44420 

.74418 

- • 44684 

.750 

.92379 

.86241 

- .44649 

.850 

.50617 

. 93293 

- .44624 

.950 

. 24706 

.96966 

-.44613 

1 .050 

.11114 

.98697 

-.44612 

1.150 

. 04738 

.99458 

-  .  4462 1 

1 .250 

.01957 

.99778 

-.44635 

1 .350 

.00797 

.99909 

- . 44652 

1 .450 

.00324 

.99963 

-.44669 

1  .550 

.00132 

.99985 

—  .44686 

1 .650 

.00055 

.99993 

- .44699 

1 .750 

.00023 

.99997 

-.44710 

1.850 

.00010 

.99999 

-.44717 

1  .950 

.00004 

.99999 

-.44721 

2.050 

.00002 

1.00000 

-.44721 

2.150 

.00001 

1.00000 

-.44717 

2.250 

.00000 

1  .00000 

-.44711 

2.350 

.00000 

1  .  00000 

-.44701 

2.450 

.00000 

1.00000 

-.44688 

ME  AN  = 

.51702  VARIANCE 

>  .04753 

NORM  = 

1 .0133 


' 
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Bj  =  -1.  0  B  =  1.  0  N  =  15 


b21 

h(b2l)/N°RM 

H(b  )/NORM 

z 

-2.450 

.00000 

.00000 

.44688 

-2.350 

.00000 

.00000 

.44701 

-2.250 

.00000 

.00000 

.4471 1 

-2. 150 

.00001 

. 00000 

.44717 

-2.050 

, 00002 

.00000 

.44721 

-1 .950 

.00004 

.  000  0 1 

.44721 

-1 . R50 

.  000  1  0 

.00001 

.44717 

-1 .750 

.00023 

.00003 

.44710 

-1 .650 

.00055 

.00007 

.44699 

-1 .550 

.00132 

.00015 

.44686 

-1 .450 

.00324 

.00037 

.44669 

-1 .350 

.00797 

.00091 

.44652 

-1 .250 

.01957 

.00222 

.44635 

-1.150 

.04738 

.00542 

.44621 

-1 .050 

.11114 

.01303 

.44612 

-.950 

.20706 

.03034 

.44613 

-.850 

.50617 

.06707 

.44624 

-.750 

. 92379 

.13759 

» 44  64  9 

-.650 

1 .44420 

.25582 

»  44684 

-.550 

1 . 85500 

.42239 

.44715 

-.450 

1 .88258 

.61227 

.44713 

-.350 

1 .46669 

.78197 

.44614 

-.250 

. 86894 

.89875 

.44307 

-.150 

. 39832 

.96074 

.43625 

-.050 

. 10787 

.98661 

.42349 

.050 

. 04750 

.99594 

.40273 

.150 

.01419 

.99879 

.37293 

.250 

.00419 

.99963 

.33471 

.35  0 

.00128 

.99988 

.29019 

.450 

.00041 

.99996 

.24212 

.550 

.00014 

.99999 

. 19309 

.650 

.00005 

.99999 

.14510 

.750 

.00002 

1.00000 

. 09949 

.850 

.00001 

1  .  00000 

. 05705 

.950 

.00000 

1  .  00000 

.01812 

1.050 

.00000 

1.00000 

-.01724 

1.150 

.00000 

1.00000 

-.04915 

MEAN  = 

-.51702  VARIANCE 

=  .04753 

NOR  M  = 

1 .01332 


' 

n 
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B2  =  0.  0  B  =  1.  0  N  =  30 

b21  h(b21)/NORM  H(b21)/NORM 


-1.700 

.00000 

.00000 

.37937 

-1.600 

.00000 

.00000 

.37519 

-1.500 

.00001 

.00000 

.37027 

-1 .400 

.00003 

.00000 

. 36445 

-1  .300 

.00009 

.00001 

. 35751 

-1.200 

.00028 

.00002 

.34920 

-1.100 

.00085 

.00008 

.33919 

-1.000 

.00255 

.00024 

.3271 1 

-.900 

.00730 

.00070 

.31248 

-.800 

.01986 

.00199 

.29478 

-.700 

.  05045 

.00534 

.27343 

-.600 

.  1  1779 

.01345 

.24786 

-.500 

.  24894 

.03129 

.21759 

-.40  0 

.46936 

. 06660 

.18232 

-.300 

.77944 

.12857 

.14214 

-.200 

1 .12828 

.22399 

.09758 

-.100 

1.41292 

.35182 

. 04969 

0.000 

1 .52373 

.50000 

0.00000 

.100 

1 .41292 

.64818 

-.04969 

.200 

1 . 1 2828 

.77601 

-.09758 

.300 

.77944 

.87143 

-.14214 

.400 

.46936 

. 93340 

-. 18232 

.500 

.24894 

.96871 

-.21759 

.600 

.  11779 

.98655 

- . 24786 

.700 

.  05045 

.99466 

-.27343 

.800 

.01986 

.99801 

-.29478 

.900 

.00730 

.99930 

-.31248 

1  .000 

.00255 

.99976 

-.3271 1 

1.100 

.00085 

.99992 

-.33919 

1 .200 

.00028 

.99998 

-.34920 

1.300 

.00009 

.99999 

-.35751 

1  .400 

.00003 

1 .00000 

-.36445 

1 .500 

.00001 

1.00000 

-.37027 

1.600 

.00000 

1 .00000 

-.37519 

1  .700 

.00000 

1.00000 

-.37937 

MEAN  = 

.00000 

VARIANCE^  .07134 

NORM  = 

1 .01403 


, 


' 
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B  =  0.  0 

b2. 

B 

c 

o 

II 

2 

h(t>2 1  )/N°RM 

H(b  )/NORM 

-1 .300 

.00000 

.00000 

-1 .200 

.00000 

.00000 

-1.100 

.00001 

.00000 

“1.000 

.00005 

.00000 

-.900 

.00022 

c  0000  1 

”.800 

.00101 

.00007 

-.700 

. 00448 

.00031 

-.600 

.01845 

.00135 

-.500 

.06819 

. 0  0532 

-.800 

.21715 

. 01866 

-.300 

.56983 

. 0564] 

“.200 

1 .17987 

. 14249 

-.100 

1 .85725 

.29516 

0.000 

2.16697 

.50000 

.100 

1 .85725 

.70484 

.200 

1  .  17987 

.85751 

.300 

.56983 

. 94359 

.400 

.21715 

.98134 

.500 

.06819 

.  99468 

.600 

.01845 

.99865 

.700 

.00448 

. 99969 

.800 

.00101 

.99993 

.900 

.00022 

. 99999 

1.000 

.00005 

1  ,  00000 

1.100 

.00001 

1.00000 

1.200 

.00000 

1.00000 

1.300 

.00000 

1.00000 

N  =  30 


z 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


MEAN=  .00000  VARIANCE^  .03571  NORM= 


1 .00837 


1 

■ 

' 
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B  =  1.  0  B2  =  0.  0  N  =  30 


b21 

h(b  )/NORM 

H(b  )/NORM 

z 

-.900 

.00000 

.00000 

0.00000 

/ 

-.800 

.00000 

.00000 

0.00000 

-.700 

.00004 

.00000 

0 .00000 

-.600 

.00046 

.00002 

0.00000 

-.500 

. 00456 

.00022 

0.00000 

-.400 

.03750 

.00195 

0.00000 

-.300 

.22913 

.01346 

0.00000 

-.200 

. 92994 

.06689 

0.00000 

-.100 

2.26299 

.22446 

0.00000 

0 .000 

3.07073 

.50000 

0.00000 

.100 

2.26299 

.77554 

0.00000 

.200 

.92994 

.9331 1 

0.00000 

.300 

.22913 

.98654 

0 .00000 

.400 

.03750 

.99805 

0.00000 

.500 

. 00456 

.99978 

0.00000 

.600 

.00046 

.99998 

0.00000 

.70  0 

.00004 

1 .00000 

0.00000 

.800 

.00000 

1.00000 

0 .00000 

.900 

.00000 

1.00000 

0.00000 

M  F  A  N  = 

.00000  VARIANCE 

=  .01756 

NOR  M  = 

1 .00634 

-  48  - 


B  =1.0  B  =1.0  N  =  30 

1  2 


b  h(bzl)/NORM  H(b21)/NORM 


-.350 

.00000 

.00000 

-.29019 

-.250 

.00002 

.00000 

-.33471 

-.150 

.00020 

.00001 

-.37293 

-.050 

.00200 

.00010 

-.40273 

.  050 

.01791 

. 00090 

-.42349 

.150 

.12262 

.00686 

-.43625 

.250 

.56241 

. 03788 

-.44307 

.350 

1 .57602 

. 14159 

-.44614 

.45  0 

2.59970 

. 35414 

-.44713 

.550 

2*56448 

.62054 

-.44715 

.65  0 

1 .59762 

. 83131 

-.44684 

.750 

.67777 

.94251 

-.44649 

.850 

.21233 

.98436 

-.44624 

.950 

. 05302 

.99643 

-.44613 

1 .050 

.01128 

.99928 

- .4461 2 

1 .150 

.00216 

.99987 

-.44621 

1  .250 

.00039 

. 99998 

- . 44635 

1 .350 

.00007 

1.00000 

-.44652 

1.450 

.00001 

1.00000 

- . 44669 

1.550 

.00000 

1 .00000 

- . 44686 

ME  AN  = 

. 50842 

VARIANCE^  .02218 

NORM  = 

1  *  0  065 


, 

-  49  - 


B  =  -1.  0  B  =1.0  N  =  30 

1  2 


b21 

h(b  )/NORM 

H  (b  )/NORM 

J. 

z 

1.550 

.00000 

.00000 

.  4468^> 

1 .450 

.00001 

.00000 

.44669 

3  .350 

.00007 

.00000 

.44652 

1.250 

.00039 

.00002 

.44635 

1 .150 

.00216 

.00013 

.44621 

1  .  050 

. 01128 

.00072 

.44612 

-.950 

.05302 

.00357 

.44613 

-.850 

.21233 

.01564 

.44624 

-.750 

.67777 

. 05749 

.44649 

-.650 

1 .59762 

. 16869 

.44684 

-.550 

2 . 56448 

.37946 

.44735 

-.450 

2.59970 

.64586 

.44713 

-.350 

1 .57602 

.85841 

.44634 

-.250 

.56241 

.96232 

.44307 

-.150 

.12262 

.99314 

.43625 

-.050 

.01791 

.99930 

.42349 

.050 

.00200 

.99990 

.40273 

.150 

.00020 

.99999 

.37293 

.250 

.00002 

1.00000 

.33471 

.350 

.00000 

1  .  00000 

.29019 

ME  AN= 

-.50842  VAR I ANCE  =  .02218 

N  0  R  M  = 

I 


1 .0065 


' 

50  - 


B1  =  °- 

°  B  2 

=  1.0 

N  =  45 

b21 

h(b  )/NORM 

H(b  )/NORM 

z 

1  *30  0 

.00000 

.00000 

.35751 

1  .?00 

.00001 

.00000 

.  34920 

1.100 

.00003 

.00000 

.33919 

1  . 00  0 

.00014 

.  00001. 

. 3P7 1 1 

-.900 

.00065 

.00004 

.31248 

-.800 

.  00P89 

. OOOPO 

.29478 

-.700 

.01159 

.00086 

.27343 

-.600 

.04100 

.  00  3P8 

.24786 

-.500 

. 1P518 

.01109 

.21759 

-.400 

. 3PP57 

.  03P56 

. 18P3P 

-  .300 

.  688 1 P 

. 08197 

.14214 

-.POO 

1  .  1960  9 

.17565 

. 09758 

-.100 

1 .67443 

. 3P  0  1  6 

. 04969 

0.000 

1 . 87461 

.50000 

0.00000 

.  1  00 

1 .67443 

.67984 

- . 04969 

.POO 

1 . 19609 

. 8P435 

- . 09758 

.300 

.6881 P 

.91803 

-. 14214 

.400 

«  3PP57 

.96744 

-.1 82 32 

.500 

.  1P518 

.98891 

-.21759 

.60  0 

.04100 

.99672 

-.24786 

.700 

.01159 

.99914 

-.27343 

.800 

.  00P89 

.99980 

-.29478 

.900 

.00065 

. 99996 

-.31248 

1  .000 

.00014 

.99999 

-.32711 

1.10  0 

.00003 

1.00000 

-.33919 

1  .  P  0  0 

.00001 

1.00000 

-.34920 

1.30  0 

.00000 

1.00000 

- . 35751 

ME  AN  = 

.  00000  VAR  I ANCE=  .04648 

norm= 

1  .  0  0  9 


* 

-  51 


B  =  0.  0 


b2  =  °.° 


N  =  45 


b  h(bzl)/NORM  H(b  )/NORM 


z 


1.000 

.00000 

.00000 

0.00000 

-.900 

.00000 

.00000 

0.00000 

-.800 

.00003 

.00000 

0.00000 

-.700 

.00028 

.00001 

0.00000 

-.600 

.00226 

.00012 

0.00000 

-.500 

.01571 

.00087 

0.00000 

-.400 

. 08761 

.00537 

0.00000 

-.300 

. 36669 

.02605 

0.00000 

-.200 

1  .  07978 

.09522 

0.00000 

-.100 

2.11 696 

. 25490 

0 .00000 

0.000 

2.66136 

.50000 

0.00000 

.100 

2.11696 

.74510 

0.00000 

.200 

1 .07978 

.90478 

0 .00000 

.300 

.36669 

.97395 

0.00000 

.400 

. 08761 

.99463 

0.00000 

.500 

.01571 

.99913 

0.00000 

.600 

.00226 

. 99988 

0.00000 

.700 

.00028 

. 99999 

0.00000 

.800 

.00003 

1  .  00000 

0  c  00000 

.900 

.00000 

1.00000 

0.00000 

1.000 

.00000 

1.00000 

0.00000 

M  E  A  N  =  .  00000  VARIANCE^  .0  2326  NORM  = 


1.0055? 


• 

-  52  - 


B  =  1.  0 


B2 


N  =  45 


b21 

b(b?  1  )/NORM 

H  (t> ^  i  )/N°RM 

z 

-.700 

.00000 

.00000 

0.00000 

-  .600 

.00001 

.00000 

0.00000 

-.500 

.00023 

.00001 

0 .00000 

-.400 

.00537 

.00021 

0 .00000 

-.30  0 

.07933 

.00343 

0.00000 

-.20  0 

.63737 

. 03360 

0.00000 

-.100 

2.39326 

. 17784 

0.00000 

0.000 

3.76885 

.50000 

0.00000 

.  100 

2.39326 

.82216 

0.00000 

.200 

.63737 

.96640 

0.00000 

.30  0 

. 07933 

.99657 

0 .00000 

.400 

.00537 

.99979 

0 .00000 

.500 

.00023 

.99999 

0.00000 

.600 

.00001 

1.00000 

0.00000 

.700 

.00000 

1.00000 

0.00000 

mfan= 

.00000  VARIANCE 

=  .01150 

NORM  = 

1 . 0042] 


I 


53  - 


B  =1.0  B  =  1.0  N  =  45 

1  2 


b  h(b  l/NORM  H(b  J/NORM  z 

21  21  21 


-*150 

.  00000 

.00000 

-.37293 

-.050 

.00007 

.00000 

-.40273 

.050 

.00187 

.00007 

-.42349 

*150 

. 03254 

.00135 

-.43625 

.250 

.31383 

.01550 

-.44307 

.350 

1.46003 

.09713 

- . 4461 4 

.450 

3.09505 

. 32830 

-.44713 

.550 

3.05653 

.65061 

-.44715 

.650 

1 .52368 

.88188 

-.44684 

.750 

.42872 

.97377 

-.44649 

.850 

.07679 

.99594 

- .44624 

.950 

.00981 

.99953 

- .44613 

1.050 

.00099 

.99995 

- .4461 2 

1.150 

.00008 

1.00000 

-.44621 

1 .250 

.00001 

1.00000 

-.44635 

1.350 

.00000 

1.00000 

-.44652 

MEAN- 

.50560 

VAR  I ANCE=  .01447 

NORM  = 

1  *0043' 


-  54  - 


B  =  -1.0  B  =  1.  0  N  =  45 


b21  h(b21)/NORM  H(bzi)/NORM  z 


1 ,350 

.00000 

.00000 

,44652 

i — ■ 

• 

PO 

'SI 

o 

.00001 

.00000 

.44635 

1 ,150 

.00008 

.00000 

.44621 

]  .  050 

, 00099 

.00005 

.44612 

-.950 

.00981 

.00047 

.44613 

-,850 

.07679 

.00406 

.44624 

-,750 

.42872 

. 02623 

.44649 

-.650 

1 ,52368 

.11812 

.44684 

-.550 

3.05653 

. 34939 

.44715 

-.450 

3.09505 

.67170 

.44713 

-  .35  0 

1 . 46003 

.90287 

.44614 

-.250 

.31383 

.98450 

.44307 

-.150 

.03254 

.99865 

.43625 

-.050 

.00187 

.99993 

.42349 

.  050 

.00007 

1.00000 

.40273 

.150  . 

.00000 

1  .  00000 

. 37293 

MEAN=  -,50560  VAR I ANCE=  ,01447  NORM= 


1  .  0  043!: 


